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Symmetric nuclear matter is studied in the self-consistent, in-medium T-matrix approach. One- 
body spectral function, optical potential, and scattering width are calculated. Properties of quasi- 
particle excitations at the Fermi surface are discussed. Dispersive self-energies are dominated by 
contributions from the 1 S'o, 35i — 3 Di, and 3 Pi partial waves. 
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INTRODUCTION 



Single-particle properties of nucleons are mod- 
ified inside the strongly interacting nuclear 
matter Q. The optical potential describes 
the average interaction of a nucleon with the 
medium. At negative energies, the single- 
particle potential reflects the modifications of 
the quasi-particles in the nuclear matter. The 
spectral function of a nucleon is a measure of 
the energy distribution of a plane wave state in 
the system, and can be observed in electron scat- 
tering experiments. 

Single-particle properties can be calculated 
from the Brueckner-Hartree-Fock (BHF) theory 

H . Hole-hole diagrams must be added to the 
lowest order results to obtain a meaningful imag- 
inary part of the self-energy. Because of the vi- 
olation of the Hugenholz-Van Hove property by 
the BHF approach, the overall normalization of 
the single-particle energies is sometimes a prob- 
lem. BHF calculations use a quasi-particle ap- 
proximation for the propagators; therefore, a fi- 
nite order diagram has kinematical limits on the 
total energy in the scattering. A finite total en- 
ergy in the interaction could introduce distor- 
tions in the far energy tails of the calculated 
spectral functions. 

Recently the in-medium T-matrix approxima- 
tion has been applied to the nuclear matter 
§ 1 H @, ! |, @- One obtains nontriv- 
ial self-energies by summing a ladder of hole- 
hole and particle-particle propagators. The self- 
energies obtained are ^-derivable (7], [IT] ; within 
this approximation, the single-particle proper- 
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ties are consistent with global thermodynamical 
properties of the system. A necessary require- 
ment for the in-medium T-matrix approach is to 
use fully self-consistent self-energies and propa- 
gators. The propagators are dressed with the 
imaginary and real parts of the self-energy, and 
the full spectral function for nucleons must be 
taken, ft requires a serious numerical effort in 
the calculations, but using a quasi-particle ap- 
proximation in the T-matrix ladder gives too 
strong pairing |5j, |s[ fl2|, too much scattering 
[|[ D, and incorrect single-particle energies 0. 



In the present work we calculate single- 
particle self-energies and spectral functions us- 
ing fully self-consistent propagators. We take 
the complete energy dependence of the self- 
energies and spectral functions into account in 
order access reliably the high energy tails of 
the spectral functions. We present results on 
the single-particle potential, scattering width 
and spectral function, obtained with a separa- 
ble Paris interaction with several partial waves 
(App. ^). In Appendix |A| we describe in details 
the numerical methods that allow us to tackle 
the problem. 



We find that normal nuclear density is at the 
limit of the region of the 3 Si — 3 D± pairing in- 
stability. Owing to the dressing of propagators 
in the gap equation the superfluid gap is very 
small, if not vanishing (depending on the details 
of the interaction in the deuteron channel). We 
leave the detailed analysis of the pairing in sym- 
metric nuclear matter with realistic interactions 
to a different work, and study in the following 
the normal state of nuclear matter. 
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II. IN-MEDIUM T-MATRIX 

In this Section we recall the basic equations 
of the approach used, more details and discus- 
sion can be found in |, § §, |l(| [0| [□], 
The T-matrix approximation sums ladder dia- 
grams with dressed particle-particle and hole- 
hole propagators for the in-medium two-particle 
propagator. The retarded T-matrix is 

(p|T(P,Q)|p')=^(p,p') 

du}idu} 2 d 3 q , . , , , 

V{p,q)A{pi,LJi)A(p2,U2) 

(q|T(P,n)|p) (1) 



32tt 5 

(i - ffa) - /M) 



ft — UJl — Ll>2 + l€ 

where pi 2 = P/2 ± q and f(uj) is the Fermi 
distribution. A partial wave expansion of the 
in-medium T-matrix is performed 

<JST) (P,n)\p)=V^ ST) {p,p) 
duJiduj2q 2 dq^(jsT) l 



(A(p 1 ,UJ 1 )A(p 2 ,UJ2))p, q 

(Wfo) -/(*»)) y , 

after angle averaging the intermediate two- 
particle propagator ((...) = / ■ ■ ■ ) 

The imaginary part of the retarded self-energy 
is obtained by closing a pair of external vertices 
in the T-matrix with a fermion propagator 

f dbJ\d 5 k 
Im£(p,u;) = / 167] ,4 A(k,uj 1 ) 

((p - k)/2|ImT(p + k, u + wi)|(p - k)/2) x 
+ + (3) 



(4) 

nu- 
real 
dis- 



where 



— 2ImE(p, u>) 



(w - p 2 /2m - RcS(p, w)) + Im£(p, w) 2 

is the self-consistent spectral function of the 
cleon and 6(w) is the Bose distribution. The 
part of the self-energy is related to ImS by a 
persion relation 



f d 



dui Im£(f>, u> 



with Yjhf{v) being the Hartree-Fock self-energy. 
Eqs. (jjj), @, (H) and (@) are to be solved itera- 
tively and at each iteration the chemical poten- 
tial p is adjusted to fulfill the condition on the 
density 



P : 



dujd 3 p 

——A{p,u})f{uj) 

l07T 4 



.16fm 



(6) 



(5) 



The T-matrix approximation sums ladder di- 
agrams contributing to the ground state energy. 
In this way it regularizes the short range core 
in the nucleon interaction similarly as the BHF 
approach. At the same time, the self-consistent 
T-matrix approximation is ^-derivable |p"lf . It 
was shown for a model interaction that if a 
fully self-consistent T-matrix calculation is per- 
formed the thermodynamic consistency relations 
are fulfilled M . In the T-matrix approximation a 
nonzero imaginary self-energy appears and leads 
to a nontrivial spectral function. It is crucial to 
keep the full off-shell dependence of the propaga- 
tors in the calculations. It is possible to treat the 
resulting energy integrals using numerical meth- 
ods described in the Appendix [A|. All the calcu- 
lations are done using the Fermi energy as the 
origin of the energy scale. The figures are also 
plotted using this convention. The calculations 
are done using a se par able parameterization of 
the Paris potential |l6| [l7|| . The Fermi energy 
obtained is — 21.4MeV and the binding energy 
— 15.1McV at normal density. This is not the 
saturation point for this interaction, but the val- 
ues quoted above are quite reasonable and give 
some confidence in the single-particle properties 
we want to study. 



III. OPTICAL POTENTIAL 

The real-part of the self-energy defines the 
single-particle pole in the propagator 

uj p = p 2 /2m + Re£(p, oj p ) . 

The free dispersion relation is modified due to 
interactions with the medium. The resulting 
effective potential Re£(p, u> p ) is attractive. It 
leads in particular to a reduction of the effective 
mass to* = We find to* ~ .85to at the 

aujp 

Fermi momentum, this value depends somewhat 
on the chosen interaction. The effective mass 
can be written |jj as the product of the A:-mass 



3 




p (MeV) 

FIG. 1: The real part of the self-energy on- 
shell (solid line), the Hartree-Fock self-energy (dot- 
ted line), and the dispersive part of the self-energy 
(dashed line). 



= 1 



m9Im£( k,uj) i 



m k /m 

mass m^/m = ^1 
mk/m ~ .62 and m^/m 



and the w- 
We find 



dlm£(fc,^) I 

t 1.37 at the Fermi 
momentum. 

The real part of the self-energy is the sum of 
the Hartree-Fock and the dispersive contribution 
(|J). The real part of the dispersive self-energy 
is negative around the Fermi energy and around 
the quasi-particle pole (Fig. |l|). The value of 
the real part of the self-energy on-shell is the 
effective potential felt by the quasi-particle in 
the medium. For positive energies of the parti- 
cle it corresponds to the optical potential. The 
depth of the single particle potential is around 
— 70MeV and it is decreasing with momentum. 
In the range of positive single particle energies 
(280MeV < p < 500MeV), the single^article 
potential can be fitted with the form 

Re£(p, lu p ) ~ iJ exp (~/3 2 fc 2 /4) 

with H = -70MeV and [3 = .58fm (dashed- 
dotted line in Fig. [j]). The range of nonlo- 
cality (3 is small reflecting the fact that the ef- 
fective mass m* is not very different from the 
free one. The depth of the potential for positive 
energies is consistent with values obtained from 
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FIG. 2: The energy dependence of the single par- 
ticle potential Re£(p, ui) at fixed momentum. 



phenomenological analyzes of the optical poten- 
tial [l9| . The single particle potential is weakly 
dependent on the temperature; at T = 4MeV 
Re£(p, ojp) is shifted up by less than lMeV from 
its value at zero temperature. 

In Fig. H is presented the energy dependence 
of the real part of the self-energy. The energy 
dependence of the single particle potential is re- 
sponsible for the time nonlocality of the opti- 
cal potential. In the whole range of frequencies 
in Fig. |^ some energy dependence of the po- 
tential can be seen. This energy dependence is 
decreasing with momentum. For energies close 
to the Fermi energy the energy dependence of 
the single-particle potential is similar to the one 
presented in Ref. [^J. The energy dependence 
determines the quasi-particle strength 



Z,, 



9ReE(p,w) 



IT 



(7) 



.73. 



At the Fermi momentum we find Z. 
On the other hand, in the range 50MeV < u) < 
200MeV Re£(p, lo) is relatively flat as function 
of energy. 



IV. SCATTERING WIDTH 

The off-shell scattering width is nonzero on 
both sides of the Fermi energy but larger at 
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FIG. 3: The single particle width — ImE(p, w) as 
function of energy for p = 0, 170, 340, and 510MeV; 
solid, dashed, dotted, and dashed-dotted lines, re- 
spectively. 



FIG. 4: The temperature dependence of the single 
particle width — ImE(p,w = 0) at the Fermi energy 
for p = 0, 170, 340, and 510MeV; solid, dashed, 
dotted, and dashed-dotted lines, respectively. 



positive energies (Fig. ||). At negative ener- 
gies the tail of the scattering width extends very 
far. It is a general feature of self-consistent cal- 
culations. Dressed propagators in the T-matrix 
ladder make the imaginary part of the T-matrix 
nonzero even at very negative energies. On the 
other hand, quasi-particle approximations have 
a kincmatical limit on the lowest energy in the 
scattering. The scattering width decreases with 
momentum. It can be understood by the fact 
that at low momenta the interaction is domi- 
nated by the strongest S waves. Other calcula- 
tions §, |, |, |, |, e], |l], H of the variational, 
BHF, Born, or T-matrix type give qualitatively 
but not quantitatively similar results. The dif- 
ferences are partly due to different short range 
properties of the interactions used. 

Consistently with general features of Fermi 
liquids at zero temperature, we find that 
ImE{p,u = 0) = (Fig. |). The same rela- 
tion, coming form the restricted phase-space, is 
fulfilled by other approximations ||, ||(], g2) . At 
higher temperatures a nonzero scattering width 
at the Fermi energy appears. As can be seen 
from Fig. |I] the scattering width at the Fermi 
energy is increasing as the square of the temper- 
ature, as expected from phase space arguments 




The large value of the scattering width above 
the Fermi energy is due to the short range part 
of the interaction potential. Splitting the total 
width into contributions from different partial 
waves shows that the 3 Si — 3 D\ partial wave is 
dominant (Fig. |J). For momenta up to twice 
the Fermi momentum and energies around the 
Fermi energy, the deuteron partial wave gives 
by far the most important contribution to the 
scattering width. It is also in this kincmatical 
region that the overall scattering width is the 
largest. We have observed a similar large con- 
tribution due to this partial wave for the Mon- 
agan Q and Yamaguchi [^5| separable poten- 
tials. Consistently, for all these interactions we 
obtain similar values for Z PF ~ .7. Large values 
of the scattering width at large energies lead to 
long tails in the spectral functions. 

The imaginary part of the self-energy at the 
quasi-particle pole is the scattering width of the 
quasi-particle. From Fig. |^ we see that around 
the Fermi momentum, narrow quasi-particle ex- 
citation exist. The scattering width behaves as 
-ImE(p,w p ) oc (p — pf) 2 , as expected from re- 
stricted phase space for scattering [^3) . The on- 
shell scattering width does not reflect the very 
large values of the off-shell scattering width (Fig. 
||) . For momenta close to the Fermi surface it is 
small; and particles with large energies have a 
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FIG. 5: The single particle width — ImE(p, uf) as 
function of energy for p = 0, 340, and 680MeV. 
The total width is plotted using the solid line. The 
dashed, dotted, and dashed-dotted lines represent 
the 3 Si — 3 Di, 1 So, and 3 Pi partial waves contribu- 
tions, respectively. 



scattering probability proportional to the total 
density and cross section. At finite temperature 
quasi-particles at the Fermi surface aquire a fi- 
nite life-time (Fig. ||). 
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FIG. 6: The single particle width — ImE(p, uj) at 
the quasiparticle pole as function of the momentum 
at T = (solid line) and T = 4MeV (dotted line). 



V. SPECTRAL FUNCTION 

The role of correlation induced by the medium 
can be judged by the modifications of the 
spectral function. Nontrivial, dispersive self- 
energies lead to broad spectral function with 
non-Lorentzian shapes and long tails. High en- 
ergy parts of the spectral functions could be re- 
vealed by electron scattering experiments. It 
is important to calculate the form of the spec- 
tral functions including contributions form short 
range correlations. In Fig. ^ are presented the 
spectral functions for three representative mo- 
menta. For zero momentum the spectral func- 
tion has a broad peak below and a background 
part above the Fermi energy. It is the reverse 
for p — 355MeV. Following the behavior of 
the scattering width, the spectral function for 
any momentum goes to zero at the Fermi en- 
ergy. For momenta close to the Fermi energy the 
spectral function has a very sharp peak at the 
single-particle energy. Its has also a significant 
background part which cannot be ignored in the 
sum rules or in the calculation of effective in- 
teractions between quasi-particles. Similarly as 
for the scattering width the self-consistent cal- 
culation gives rise to a long tail in the spectral 
function at negative energies. 
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FIG. 7: The spectral function A(p,co) as function 
of energy for p — 0, 255, and 350MeV; solid, dashed, 
and dotted lines, respectively. 



FIG. 8: The momentum distribution of nucleons 
at zero temperature (solid line) and at T = 4MeV 
(dotted line). 



The spectral function can be used to obtain 
the momentum distribution 

n(p) = J^A(p,u > )f(u). (8) 

Short range correlations modifying the spectral 
function are reflected in the nucleon momentum 
distribution. The free Fermi distribution is de- 
pleted below pp and a high momentum tail in 
n(p) is formed. At the Fermi momentum the dis- 
continuity in the occupancy is reduced from 1 to 
Z PF (Q). The Fermi momentum in the interact- 
ing system should be the same as in the free one 
|fll| , |2q| . This relation is approximately fulfilled 
in the present calculation, but not as well as in 
our previous work using only S wave interactions 
0. A possible source of the discrepancy could 
be the use of the partial-wave expansion which 
spoils the <I>-derivability of the self-energies. At 
finite temperature the Fermi surface is washed 
out, but the depletion of the Fermi sea and the 
high momentum tail are the same. It means that 
the short range correlations stay essentially the 
same at T = 4MeV. The same can also be seen 
by inspecting the self-energies. At finite tem- 
perature the scattering width is increased only 
in the vicinity of the Fermi energy. 

Long tails in the spectral function have impli- 
cations for saturation properties of nuclear mat- 



ter. Increased kinetic energy due to the high mo- 
mentum tail in n(j>) is compensated by increased 
removal energies due to the negative energy tail 
of the spectral function. A detailed discussion at 
different densities will be presented elsewhere. 



VI. SUMMARY 

We study the single-particle properties of nu- 
cleons in nuclear matter using a conserving, in- 
medium T-matrix approximation. The calcu- 
lations are done in a self-consistent way with 
dressed propagators in the ladder diagrams. To 
our knowledge, it is the first such calculation 
in the literature using realistic interactions and 
several partial waves. The full discretization of 
the spectral function and self-energies allows to 
discuss the details of their energy dependence. 
We find that the basic features of a consistent 
approximation to a Fermi liquid system are ful- 
filled. The scattering width is zero at the Fermi 
energy, and the quasi-particles at the Fermi sur- 
face have infinite lifetime. The momentum oc- 
cupancy has a discontinuity of Z PF at the Fermi 
momentum. The off-shell scattering width is 
very large at energies ~ 300 — 400MeV and mo- 
menta below 400MeV. In this region the main 
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contribution comes from the 3 5i — 3 D\ par- 
tial wave. The 3 Pi and 1 S partial waves give 
also important contributions to the on-shell self- 
energies. The scattering width on-shell is not 
very large, maximally 26MeV. Finite tempera- 
ture effects are concentrated around the Fermi 
surface. At finite temperature, the scattering 
width gets finite around the Fermi energy with- 
out modifying the real single-particle potential 
and the short range correlations. The quasi- 



particle renormalization factor is Z„ 



.73 and 



is largely independent on the interaction used. 
The effective mass is m* ~ .85m. 
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APPENDIX A: NUMERICAL METHODS 

Calculations using self-consistent spectral 
functions require the evaluation of energy inte- 
grals in Eqs. (Q) and (|^). This is the main nu- 
merical difficulty for any self-consistent approxi- 
mation; only in the last years self-consistent ap- 
proaches have been applied to the study of high 
T c superconductors |27l |2^, ^9|, |3^] and nuclear 
matter [Q, ||, @, ||, §| ■ Some calculations are per- 
formed in the imaginary time formalism |2^, p£| 
which requires a numerical procedure for the an- 
alytical continuation to calculate the spectral 
function. A simpler way is to use the real-time 
formalism operating with the retarded T-matrix, 
the retarded self-energy, and the spectral func- 
tion n 
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The real-time formalism was also 
used for other calculations of the nuclear matter 
H [| performed at zero temperatures. To deal 
with the off-shell propagation a numerical pa- 
rameterization of the energy dependence of the 
spectral function A(p, uj) by a set of Gaussians 
has been used |). Alternatively, the spectral 
function can be represented as a sum of S func- 
tions |], 0- The above two methods can be 
easily applied at zero temperature, where a nar- 
row Gausian or one of the 5 functions describes 
the quasi-particle peak for momenta close to the 
Fermi surface. 

The alternative approach, here employed, uses 
a direct discretization of the spectral function 



and self-energies as functions of the energy. If 
the discretization in uj is uniform the energy in- 
tegrals in Eqs. (Q) and (|^) can be performed 
efficiently by Fourier transforms. This method 
can be used directly at finite temperature [js], |2{| , 
but not for low temperatures where the quasi- 
particles are very narrow around the Fermi sur- 
face (Fig. |^). At zero temperature we have 

A{p F ,uj) = 2irZ PF 6(uj-E F ) + B{p,uj) (Al) 

where the background part B is smooth and can 
be discretized. Generally, for any momentum 
close to the Fermi momentum the spectral func- 
tion is rapidly changing in the vicinity of the 
quasi-particle peak; therefore, it is has to be split 
into a quasi-particle peak and a smooth back- 
ground. The sharp peak is approximated by a 5 
function and the smooth part is discretized. We 
have found numerically that this separation in 
the spectral function is needed if the width of 
the quasi-particle peak — 2Im£(p, uj p ) is smaller 
than ~ 3 times the spacing in the energy dis- 
cretization Auj. We use the following represen- 
tation of the spectral function 

A(p, u) = 2ttZ p S(lo - uj p ) + B(p, uj) . (A2) 

The background part is defined as 



B(p,w) = 



-2ImS(p,ui) 

for (uj — uj p ) 2 + Im£(p, uj) 2 > r) 

-2Im£(p, uj)/t] ; 
for {uj - ujp) 2 + JmE(p, uj) 2 < rj 

(A3) 



The parameter rj is set to cut off the rapidly 
changing part of the spectral function. We use 
r\ ~ 4 Aw 2 . The strength of the quasi-particle 
component Z p in Eq. (A2) is adjusted to con- 
serve the sum rule J ^A(p,uj) = 1. The weight 
of the singular part Z p is not the same as the 



renormalization factor Z p (Eq. ^|). In Fig. 
is shown an example of the separation of the 
background and singular parts of the spectral 
function A(p = 255MeV, uj). 

For a separable interaction V a (p,p ) = 
T^i j ^f.j9i(p)9f(P ) the T a -matrix for the par- 
tial wave a = (JST) is {p\T a (uj, P)\p) = 
Y,i. 3 T l a j(u,P)g?(p)gf(p')- Where the matrix 
T%'= (Xfj-Jtj)' 1 with the coefficients given 
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by the integrals 




co (MeV) 



FIG. 9: The spectral function A(p, to) as function of 
the energy (solid line) for p — 240MeV and the corre- 
sponding smooth background part B(p,ui) (dashed 
line). The position of the S function representing the 

quasi-particle peak at uj = -7.7MeV is indicated by Substituting (A2) for the spectral functions we 
the dotted line. have 



(i - 2/M) B{ P2 ,n u) 

/d?k 
— g °(k)g°(k)B( Pl ,n- lo P2 )z P2 (i - /(n -lo P2 )- /kj) 



/ prf(W) ^ (1 " /fal " + - M • (A5) 



The en ergy integral in the first term of 
Eq. ( A5) is a convolution of the functions 
B(pi, u) (1 — 2f(u>)) and B(p 2 ,uj)- It can 
be evaluated by Fourier transformations using 
FFT transform algorithms Moreover, the 

Fourier transform and its inverse can be per- 
formed outside of the mom entum integral |3^ j. 
The second term of (A.5) is a standard two- 



dimensional integral witho ut s ingularities. The 
integral in the last term of ( |A5| ) is of the type of- 
ten occurring in quasi-particle calculations, such 
as the G-matrix or the quasi-particle T-matrix 
approximations. Angle averaging separately the 
numerator and the denominator under the inte- 
gral we have 
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Im 



< ^Pi^P2 (! - /( w Pi) - /( W P 2 )) >P,k 



k 2 dk 



fl— < LO 
-'Pi "^P2 



Pi 



U P2 



>p,k 



kl < Z m Z Vo (1 - f(uJ Pl ) - /(Wpj,)) >p,fe D 



27r a < + w P2 >p, fc /dk\ k=ka 
I 



(A6) 



where /co is the solution of 

SI =< uj v , + w m 



>P,k Q 



The angle average in the denominator < ui Pl + 
lu P2 > p t k is a function of the total and relative 
momentum (P and k). It can be represented 
using a one-dimensional function 



< Upi + c^ P2 > P>fe = — (P(P/2 + fc) - F(P/2 - k)) 



(A7) 



with 



F(.t) = / pdpuj p 



This effectively one-dimensional parameteriza- 
tion allows to perform the integral corresponding 
to the quasi-particle part of the spectral function 
very efficiently without using any parabolic ap- 
proximation for ujp. The same method can also 
be used in G-matrix calculations. The real part 
of the integrals is obtained using the disper- 
sion relation 



ReJ2(n,P) 



ft' ImJ°(»>) 

~ n' - n 



The calculation of the energy integral in the 
equation for the self-energy (|3|) proceeds very 
similarly. We have 

ImS(p,w)= J B(k,uji) 

<(p - k)/2|ImT(p + k, w + wi)|(p - k)/2) A 
(/(wi) + 6(« + «i)) 
d 3 k 
8^ 

((p - k)/2|ImT(p + k, u + w fc )|(p - k)/2)^ 
(/K) + 6(w + w fc )) ■ (A8) 



The energy integral in the first term is again 
of the form of a convolution and is calculated 
using Fourier transforms. The integral in the 
second term is a standard two-dimensional inte- 
gral. Unlike in the G-matrix approximation, the 
calculation of the self-energy requires the knowl- 
edge of the full off-shell T-matrix. Restricting 
oneself to contributions from on-shell T-matrix 
leads to erroneous results |3^] . 

In the iteration of the self-consistent set of 
equations (|l|, [| |J, H an d ||) it is advantageous to 
parameterize the energy dependence of the off- 
shell quantities with respect to the Fermi energy. 
The iterations are much faster, since a modifica- 
tion of the chemical potential between iterations 
does not change the most important features of 
the u dependence of E(p, w) and T(Q,P), e.g. 
ImE(p, lu = 0) = and the pairing singularity 
for T(fi = 0, P = 0). The absolute energy scale 
is recovered only in the calculation of the to- 
tal density (@) and final observables. This way 
of proceeding with the iterations can make use 
of self-energies calculated for other densities or 
temperatures to start a new iteration. 



APPENDIX B: N-N INTERACTION AND 
NUMERICAL PARAMETERS 

We use a separable parameterization of the 
Paris potential fl6, 17|. It contains all the J = 0, 
J = 1, J = 2, and the 3 D^ — 3 G3 partial waves. 
For the most important So and 3 Si — 3 D\ par- 
tial waves we choose the rank 3 and 4 respec- 
tively. We have observed that the use of the 
separable parameterization of the Paris poten- 
tial for the 3 Pi partial wave leads to very high 
values of the effective mass m* ~ .95m. Since 
the 3 Pi phase shifts are not well reproduced by 
the parameterization of Ref. |lq] , we choose the 
Mongan / parameterization |24|] for this partial 
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wave. It is important to check that no unphysi- 
cal bound states occur in the off-shell T-matrix 
for the given choice of N — N potentials. The 
contribution from such unphysical bound states, 
even far from the considered on-shell energies, 
would spoil the calculated self-energies. 

The use of the Fourier transform algorithm 
for the calculation of energy integrals requires a 
fixed spacing in the energy grid. It means that 
we have to set finite ranges for the kinematical 
variables in the model. The single particle mo- 
menta are in the range p < 1700MeV. The to- 
tal momentum of a nucleon pair is limited by 



3400MeV. The energy dependence of the T- 
matrix T(fi, P) is calculated for -2500MeV + 
P 2 /Am < Q < 3500MeV. The energy range 
for the single particle self-energy is taken from 
-1700MeV to 3500MeV. The energy dependent 
functions are discretized with a grid spacing of 
1.28MeV. 

To obtain results with stability between iter- 
ations better than 1% we need typically 7 iter- 
ations at zero temperature. The convergence of 
the iterations is faster and more stable at finite 
temperature. 
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